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ff^ • Abstract. Thermodynamic Bethe ansatz equations are coupled non-linear 

1^^ ' integral equations which appear frequently when solving integrable models. 

I ^ Those associated with models with N=2 supersymmetry can be related to 

, differential equations, among them Painleve III and the Toda hierarchy. In 

C ' the simplest such case the massless limit of these non-linear integral equations 

^2 ' can be solved in terms of the Airy function. This is the only known closed-form 

^~i ^ solution of thermodynamic Bethe ansatz equations, outside of free or classical 

|/^ , models. This turns out to give the spectral determinant of the Schrodinger 

^_^ ' equation in a linear potential. 

1—i '. 

^ , A great deal of interesting mathematical physics has arisen from the study of 

,_^ ■ integrable models of statistical mechanics and field theory. One interesting area is 

y—i , known as the thermodynamic Bethe ansatz (TBA), which has proven a useful tool 

^»0 ' for computing the free energy of an integrable 1-1-1 dimensional system |l| . One ends 

^^ I up with a set of coupled non- linear integral equations, the "TBA equations". One 

zi^ ' completely- unexpected result was a correspondence between a limit of these integral 

■"^ , equations and some very well-studied non-linear differential equations, namely the 

r^ • Toda hierarchy 0. The purpose of this paper is to extend these results further, 

and show that in at least one case there is a closed-form but non-trivial solution of 

the integral equations. Not only is it interesting that such complicated equations 

(-* \ have a simple solution in terms of the Airy function, but proving it requires some 

?• ' utilizing some very intricate results involving the Painleve III differential equation 

.^H ! HI ^, §]■ Moreover, it turns out to be related to the spectral determinant of the 

J^ ' Schrodinger equation in a linear potential [g], |6|, |7|. 

H I The TBA integral equations are generically of the form 

(1) eaie) = ma cosh^ -J2[ ^^ab{9 - 9') ln(l + e^^-^^C'')). 

Physically, Tea{0) is the energy for creating a particle of type a and rapidity in a 
thermal bath at temperature T. The nia are the particle masses over temperature, 
while the fia are their chemical potentials over temperature. The kernels (/>afc are 
a result of the interactions between particles. This and all unlabelled integrals in 
this paper run from — oo to cx). The free energy per unit length is 

(2) F = -T'^^ I — ma cosh 6* ln(l + e^^-^^-^^)). 

Here we study TBA equations where the underlying physical system has N=2 
supersymmetry. The amazing thing is that solutions of a particular limit of such 
TBA equations are simply related to solutions of a non-linear differential equation 
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Oh! 
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I2I. Particles in an N—2 theory all have a charge, /a, known as their fermion num- 
ber. When the chemical potentials are /Iq = ^'Kfa^ a consequence of supersymmetry 
is that the Eq in (|]) are all constants and the free energy is F = 0. This is known as 
the Witten index (the usual integer contributions to the index do not contribute to 
the free energy per unit length) Q . The result of g] is that for chemical potentials 
Ma — ii^^ ~ h)fa, one can derive a differential equation for the order h piece of the 
free energy. The simplest cases give the Painleve III differential equation and the 
Toda hierarchy. 

The TBA equations for the case at hand were derived in M. They have (j)i2 = 
013 = l/cosh0 with (j)ab = ipba and other (ji^b = 0; while ?Ti2 = "^3 = 0, and 
-/3 = 1. For small positive h, the functions A and B are defined by 
- In/i, and €2{0) = £3(6*) = —hB{9). The order h TBA equations are 



/i=0 

6l(0) = 


h = 

A{e) 


(3) 




(4) 





Bi,) ^ I 



^^' ^ r~A(e') 



27r cosh(6l - 6') 

where 2u{9) = mi cosh 6* here. In |g], a physics proof was given that the resulting 
free energy is simply related to a solution of the Painleve III differential equation 
with variable mi. This is a physics proof because one method of computation 
gives integral equations, the other the differential equations. This result was ex- 
tended considerably in [g|. Subsequently, the equivalence was proven directly and 
rigorously in jj] . 

A particularly interesting situation is the "massless" limit, where m,i is very 
small. Then 2u{9) = e^, because mi can be removed by redefining by a shift. 
The result of this paper is A and B in (^J4|) can be found in closed form in this 
massless limit. 

Result: 

When u{e) = 6^2, the solution of (|[|) is 

(5) e-^(^) = -27r^(A*(z))2 

dz 

(6) B{9) = -27r-^ Ui(ze"/3)^j(;je-"/3) 

dz L 

where z = (3e®/4)^/'^ and Ai{z) is the Airy function. 

To check the normalization, note that (|,|) imply that e~^(^) -^ 2/\/3 and B{e) -^ 
l/vS as 9 —> —00, in agreement with the limits of the appropriate Airy functions 
as z — > 0. Ai{z) is a solution of the differential equation w" = zw, while e~"^*-^) 
and B{9) solve w'" — Azw' ~ 6w. 

As far as I can tell, it is difficult to prove the result by direct substitution into the 
integrals, but requires utilizing some additional structure. First, define the integral 
operator K which maps functions to functions with the kernel 

e" + e" 
where 
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Solutions of Painleve III can be expressed in terms of this operator K [|l^ || || 
when 2m — mi coshfl. The functions Z+ and Z_ are defined as 

(7) Z+ = 6-29/3(7 + ^)-i^ Z_=e-'^/\l~K)-'^E 

These Z± are simply related to the Q^ P — (/± K)^^E used in Q. It was proven 
in |, |, ll that functions A and B solve (|[|) if 

(8) e-^C') = ATrZ+{e)Z_ie) 

(9) B{9) ^ 47^e''^/3e"(^+''^/2)+«(e-^V2)^_^(g)_^^/2)Z_(6l + ^7^/2)-i 

The functions (|,|) obey (|,|) for any entire u{e) obeying u{9) = u{9 + 2m) §, ||, 
not only the special case u{9) — e^/2 analyzed in this paper. 

With the similarity between (Hj) and (||^), proving the result is equivalent to 
proving that 

(10) Z+{e) = Ai{z) Z_{0) = -Ai'{z) 

when u{6) = e^/2. The expression (0) for B{9) follows by using a few Airy-function 
identities. The expressions for Z± in (nfl) can be proven by directly evaluating the 
integral in their definition. Specifically, ([^) requires that 

(11) E{9) = e2«/3z+(0) + 2E{e) / d9'Z+{e')e-'' /^ _^ 



Defining x = e^/2 and using the expression of an Airy function in terms of a Bessel 
function, the integral is proportional to 



v2x 

This can be looked up in [0|, or evaluated by using Mellin transforms, under which 
the Bessel function K^ has nice properties. One indeed finds that (O) and the 
analogous relation for Z- are true when Z± are given by (M). Since the solution 
with the appropriate analyticity properties (no zeroes and bounded in the strip 
\Im9\ < 7r/2) is unique [Q, the result is proven. 

The functions Z+{9) and Z-{9) are interesting in their own right. They arose as 
a sort of partition function of the boundary-sine Gordon problem |l3] at coupling 
(7=2/3, and equivalently as the continuum version of the Baxter Q-operator |l3] in 
studies of conformal field theory [|l^ . The Q operator gives the generating function 
of the conserved quantities (local and non-local) of the theory. Up to normalization, 
the functions Z± are Zbsg{z, ±1/3) in H, and (Qiz)) at p = ±1/6 in [|4|. Thus 
the result here provides the only case (other than where the system is free or 
classical) where these quantities can be computed explicitly. In fact, it provides 
strong evidence that the results of [|l^, |4| are applicable in the repulsive regime 
g > 1/2 as conjectured. For example, it shows that the zeroes of the eigenvalues of 
the Q-operator obey the pattern conjectured in [h4|. The "quantum Wronskian" 
of [Q becomes equivalent to the Proposition of^, and both are easily verified 
by using the (ordinary) Wronskian of the Airy functions. Also, by using the series 
expansion of |1^ , it gives a strange sequence of identities of sums of products of 
gamma functions. 



Recently, the work of |14 also arose in a completely new context. Namely, it 



was observed in M that the spectral determinant for the Schrodinger equation in a 
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potential |a;|" for any a is given precisely by the vacuum eigenvalue of this Baxter Q- 
operator. A non-zero angular momentum can be added, and this correspondence 
still holds 0. This is because the quantum Wronskian is equivalent to a set of 
functional relations for the spectral determinant derived in ^, . The result above 
then means that Airy function is the spectral determinant for the Schrodinger 
equation in a linear potential, a result shown directly in [p|. It would be most 
interesting to understand how to extend this result. 

Since this system discussed in this paper provides the simplest example of the 
differential equation/TBA correspondence , it seems likely that there is a simple 
solution of the massless limit of any TBA equations of this type. What is not yet 
known outside of this case is the detailed analysis of the differential equation of iQ , 
which was vital to the results of [El . Given how beautifully the Airy function solved 
the problem here, it would be quite interesting to see how this result is generalized. 

This work was mostly done in spring 1996, while I was a postdoc at the Uni- 
versity of Southern California. I thank Hubert Saleur and the Packard Foundation 
for support then. It was presented in 1997 at the Supersymmetry and Integrable 
Models Workshop, University of Illinois-Chicago, and the International Workshop 
on Statistical Mechanics and Integrable Systems in Coolongatta. I thank the or- 
ganizers for their hospitality. I thank P. Dorey for conversations, and S. Lukyanov 
for urging me to write this up despite the late date. Currently, my work is sup- 
ported by a DOE OJI Award, a Sloan Foundation Fellowship, and by NSF grant 
DMR-9802813. 
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